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Abstract 

Essential properties of multilayer graphenes are diversified by the number of layers 
and the stacking configurations. For an A-layer system, Landau levels are divided into 
N groups, with each identified by a dominant sublattice associated with the stacking 
configuration. We focus on the main characteristics of Landau levels, including the 
degeneracy, wave functions, quantum numbers, onset energies, field-dependent energy 
spectra, semiconductor-metal transitions, and crossing patterns, which are reflected 
in the magneto-optical spectroscopy, scanning tunneling spectroscopy, and quantum 
transport experiments. The Landau levels in AA-stacked graphene are responsible for 
multiple Dirac cones, while in AB-stacked graphene the Dirac properties depend on 
the number of graphene layers, and in ABC-stacked graphene the low-lying levels are 
related to surface states. The Landau-level mixing leads to anticrossings patterns in 
energy spectra, which are seen for intergroup Landau levels in AB-stacked graphene, 
while in particular, a formation of both intergroup and intragroup anticrossings is ob¬ 
served in ABC-stacked graphene. The aforementioned magneto-electronic properties 
lead to diverse optical spectra, plasma spectra, and transport properties when the 
stacking order and the number of layers are varied. The calculations are in agreement 
with optical and transport experiments, and novel features that have not yet been 
verified experimentally are presented. 
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1 Introduction 


Graphene is a one-atom-thick material made up of carbon atoms arranged in a hexagonal 
lattice structure HM. The particular electronic properties and the high carrier mobility 
(up to 15,000 cm^V“^s“^) (SHH] make this new material a promising candidate for next- 
generation nano-devices [TOHTS]. Graphene is the building block of other carbon allotropes, 
geometrically in terms of curling, cutting, folding and stacking— such as the synthesized 
zero-dimensional (OD) fullerenes [TB], ID graphene nanoribbons [TTlfTB] and carbon nan¬ 
otubes [20], 2D layered graphenes BE] and 3D bulk graphite [21]. These low-dimensional 
systems have been a focus of intense interest in the areas of electronic properties [221 - I52] . 
optical properties |i6] - H8l[55|l6Tl [65 ]l86flll3] . Goulomb excitations [12211133] and quantum 
Hall transport properties |^ [30|[TMlll46] . It has been demonstrated in theoretical and 
experimental research studies that the essential properties are dominated by the inter¬ 
play between the external helds and the specihc geometric symmetries, and related to 
the geometric boundary conditions [221 - IBT] . The results are rather different from those in 
conventional metals and semiconductors [14711148] . To illustrate the geometry-dependent 
response under magnetic quantization, we thoroughly review multilayer graphenes with 
different stacking conhgurations. 

Planar monolayer graphene is a starting point toward the broad investigation of other 
carbon allotropes with diverse geometric symmetries. It was a theoretically studied mate¬ 
rial until in 2004 a group led by Andre Geim and Kostya Novoselov extracted monolayer 
graphenes from bulk graphite by using the micro-mechanical exfoliation technique [1] . Their 
success in isolating graphene sheets immediately triggered much activity in obtaining high- 
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quality large graphene sheets. As a result of the efforts of many groups over the past 
decade, the improved quality is adequate for research experiments and industry applica¬ 
tions [^ ll49f[T^I203] ■ The electronic structures of graphenes were hrst derived by Wallace 
from the simple nearest-neighboring tight-binding model jl53j . The calculated results show 
that monolayer graphene is a gapless semiconductor with a vanishing density of states at 
the Fermi level {Ep = 0). The low-lying energy dispersions linearly depend on the wave 
vector, with crossings of conduction and valence bands at the K and K' points (the hexag¬ 
onal corners of the hrst Brillouin zone). This leads to a conical energy spectrum around 
the corners K and K' of the Brillouin zone, resembling the Dirac fermion spectrum of 
relativistic quantum particles, with an effective speed of light ~ c/300, where c is the light 
speed. As a result, the quasi-particles in graphene can be described by a Dirac-like Hamil¬ 
tonian HEllET]. The conical structures are known as Dirac cones, and the similarly band 
crossing points, K and K\ are known as Dirac points. 

Graphene sheets are held together by weak Van der Waals interactions. The stacked 
graphenes have been experimentally characterized by infrared spectroscopy [981 - 1101] , trans¬ 
mission electron microscopy [15411158] . scanning tunneling microscopy (STM) and spec¬ 
troscopy (STS) [159fll71| . Raman spectroscopy [175fll81| . and angle-resolved photoelectron 
spectroscopy (ARPES) [18211188| . The sacking orders are arranged in specihed sequences, 
essentially AA AB [157[ll59[ll001[T75lfT79] and ABC stacking [THTlfTonifTTmTTO] . 

A haphazard stacking of graphene sheets is called turbostratic structure [18011181] . AB and 
ABC stacking conhgurations are common stacking orders of natural graphite [18911190] . 
while AA stacking conhguration has been found only in intercalated graphite [191] . It was 
reported that natural graphite typically contains 80 % of volume fraction of AB conhgura- 
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tion, 14 % of ABC configuration, and the rest part is about the other turbostratic struc¬ 


tures |189lll90j . Many works on graphene synthesis have been devoted to obtain few-layer 
graphenes by exfoliation of highly orientated pyrolytic graphite (mechanical and chemical 
processes) [nil92pi93j . chemical and electrochemical reduction of graphene oxide |194fll98j . 
and arc discharge |199H201] . The process of mechanical exfoliation is low-yield, size limited 
and cannot produce large scale graphenes. However, chemical methods have the advantage 
of yielding large-scale graphenes, but the drawback is the inferior quality of the exfoliated 
graphene sheets which contain many defects during the oxidation and reduction processes. 
Recently, chemical vapor deposition method promises as a potential means for large areas 
of high quality graphenes onto the transition-metal substrates (e.g., Ni, Pd, Pt, Ir, and 
Cu) based on the saturation of the hydrocarbon gas and the substrate [154111581120211207] . 
Previous works conhrmed all the three highly symmetric stacking conhgurations, AA, AB 
and ABC, in the CVD-grown few-layer graphene [15411158] . Although AA conhguration 
is theoretically predicted with higher total energy [IT], the layer structure was grown on 
on the (111) surface of diamond [154] and on the SiC (OOOi) substrate [155] . Moreover 
this structure was frequently observed in bilayer graphenes [156] . On the other hand, AB 
structure is mixed with ABC one [1581ll60lll7711T79] since the tiny difference between the 
total energies of the two conhgurations makes the difficulty of selectively producing AB or 
ABC stacking [4711157 . 

AA-stacked graphenes have the simplest geometric structure, where the x- and y- 
coordinates of carbon atoms on different layers are vertically aligned [33lll541[155] . Their 
band structure preserves the conical structure due to the vertically projected geometric 
symmetry; several Dirac cones are exactly located at the K {K') point with energy separa- 
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tions in the order of the vertical interlayer hopping energy [321ISSIISS] • In the AB stacking 
conhguration, the second graphene sheet is shifted by one bond length along the armchair 
direction with respect to the first layer [39] . If the third layer is further shifted in the same 
way, the arrangement sequence is known as an ABC stacking [391 • AB- and ABC-stacked 
graphenes show more distinct band structures than AA ones as a result of the signihcant 
hybridization of energy bands under the interlayer atomic interactions [351133113911301133] . 
AB-stacked graphenes consist of several bilayer-like parabolic bands, while it owns a partic¬ 
ular pair of monolayer-like linear bands if the total number of layers is odd [?.I35 ]I36]I38[ I39]. 
In ABC systems, there is always one pair of weakly dispersive bands with surface-localized 
states near the Fermi level [3911331133935] . and peculiar sombrero-shaped dispersion bands 
near the energy of the vertical nearest-neighboring interlayer interaction [39932] . Several 
Dirac points are located at low energies for both AB- and ABC-stacked graphenes with 
their characteristics related to the number of layers and the stacking. The above-mentioned 
stacking-dependent electronic structures have been experimentally verified by infrared op¬ 
tical spectroscopy [989102] . and angle-resolved photoemission spectroscopy [1829188] 

The electronic properties under the influence of external helds is focal point of some 
research [331133113711391139115911531 IHm 1939m ^ In the presence of a uniform 

perpendicular electric held, the Dirac cones in AA-stacked graphenes are preserved [331133] 
l65] . while only a rigid shift from Ep = 0 takes place. However, the electric held can induce 
a band gap in few-layer AB- and ABC-stacked graphenes [39953] . Optical spectroscopy 
[53957] and transport experiments [58960] conhrmed that the gap can be continuously 
tuned up to hundreds of meV. An tunable band gap makes them interesting for wide range 
of applications of electronic and photonic devices, such as electronic transistors [IHllII] and 
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photodetectors [IMS]. 

Electronic states in a uniform perpendicular magnetic field B = Bqz are evolved to 
quantized Landau levels (LLs) In density of states (DOS), the peculiar sequence 

of peaks is responsible for the quantization of massless Dirac carriers USSHIZD]; this se¬ 
quence is dependent on \/nBo, where n is the quantum number and Bq is the magnetic 
held strength [11 I134[[T^ . On the other hand, the massive fermions in bilayer graphene 
give rise to quantized LL energies [671ll39[[T7r| . The held-independent LL at 

Ep = 0 for n = 0 conhrms that the Dirac points are hole-electron states in charge neutral 
graphene (available to either holes or electrons). The LL spectra of multilayer graphenes 
are particularly depend on the layer number and stacking conhguration [63H82]. Based on 
the zero-held band structures, the LLs can be classihed into several groups, with each be¬ 
ing responsible for the magnetic quantization of the respective subband. The nature of the 
Dirac points is the origin of the LLs with constant energy in monolayer graphene [LllH ll35j 
and with held-insensitive energies in few-layer graphenes [651176] . It should be noted that 
such LLs in AA-stacked graphenes are separated from each other by tens of meV [331166] . 
while such LLs in AB- and ABC-stacked graphenes are conhned in a narrow energy range 
of 10 meV around Ep = 0 [671I691I7311751I7611H2]; their energy diherences can be verihed 
by quantum transport experiments |138[I141]I143[I144] . Furthermore, one LL is coupled to 
another with a specihc quantum number relationship, leading to the anticrossing patterns 
in LL spectra, which has been investigated in few-layer AB- and ABC-stacked graphenes 
and in the Bernal-stacked bulk graphite [6^1781181118^1208] . This special phenomenon can 
be explained by a perturbation of the non-vertical interlayer atomic interactions [208] . 
Previous works based on the comprehensive tight-binding model [75l|82] showed that the 
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Landau quantization of the sombrero-shaped subbands in ABC-stacked graphenes leads 
to the peculiar LL evolution under a variation of the magnetic held that reveals inversely 
held-dependent energies and a complex pattern of intragroup LL anticrossings, especially 
in the region of weak helds. 

Scanning tunneling spectroscopy (STS) has been used to study the Dirac-fermion prop¬ 
erties of graphenes through the measurements of the DOS |159H165] . In monolayer graphene, 
the DOS is linearly dependent on the energy and vanishes at the Dirac point, a behavior 
identical to that of the massless Dirac fermions |161[I162] . It has been verihed that massless 
and massive Dirac spectra can appear in AA- and AB-stacking bilayer graphenes, respec¬ 
tively |162] . Furthermore, the coexistence of both spectra can be observed in symmetry- 
broken bilayer graphenes |162[|163j . In a magnetic held, the LL energies can be directly mea¬ 
sured through the energies of the sharp peaks in the diherential conductance map of STS 
{dl/dV versus V) [16611171] . The orbital Landau quantizations of massless Dirac fermions 
in monolayer [16611170] and massive Dirac fermions in AB-stacked bilayer [171] are, respec¬ 
tively, verihed by the sequence of peaks with square-root and linear dependencies on Bq. 
In terms of the sequence of peaks, the indexes n = 0 and n = 0,1 correspond to the mono- 
layer and bilayer Dirac points, respectively. On the other hand, the spectroscopic-imaging 
STM has been used to obtain the local DOS (LDOS) of LLs, such as in graphene [172] . 2D 
electron gas [173] and topological insulator [174] . revealing the spatial distributions of the 
wave functions which consist of n-dependent nodal structures [17311174] . Owing to the LL 
anticrossings and the symmetry broken down in multilayer graphenes [38ll69ll75ll771182[|136] . 
it is interesting to measure the behavior of LLs and wave functions. 

Optical spectroscopy can be utilized to identify the feature-rich electronic properties 
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of graphenes [Ml - IFn[^5f4112j . The optical response is determined by vertical transitions 
from the occupied states to the unoccupied states. Experimental [9711102] and theoreti¬ 
cal [88H92] studies are consistent in the spectral absorption features of multilayer graphenes 
without external fields, such as the prominent structures, and the number, intensity and 
frequency of absorption peaks, which depend on the layer number and stacking conhg- 
uration. Furthermore, the theoretically predicted gap opening of biased AB-stacked bi¬ 
layer and ABC-stacked trilayer graphenes has been verified through the observation of 
new transition channels within the energy gap [SlHFT] . The magnitude of the energy 
gap is compatible to that calculated within the tight-binding model over a wide spectral 
range [l9H53l[92]. The magneto-optical absorption frequencies, identified as monolayer and 
bilayer graphene features, show a square-root dependence [10411108] and a linear depen¬ 
dence ptiUTl fTHmiTT^ on i?o, respectively. The square-root dependent frequencies are also 
expected for the LL transitions of Dirac quasi-particles in the graphene-like 2D materials, 
e.g., M 0 S 2 [11411115] and silicene [11611117] and the topological insulator [11811121] Accord¬ 
ing to the Fermi golden rule, the inter-LL transitions from the well-behaved LLs obey the 
specific selection rule An = ±1 [9411971 [T06llll411TT6lll201[T^ . With an increase of the 
layer number, the magneto-electronic properties are enriched and diversified by distinct 
stackings. A further exploration on the multilayer spectra needs to be undertaken in order 
to comprehend the magneto-optical response to the special LL spectra. 

In electronic transport measurements, the peculiar Landau quantization of massless 
Dirac fermions results in equidistant plateaus of the Hall conductivity of axy = ±4(n -1 
\/2)e^/h, where e ^/h is the conductance quantum and n is an integer number [5E1I134H138] . 
The quantum Hall effect (QHE) can be observed even at room temperature due to the high 
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mobility of electrons in graphene jl36j . Unlike the sitnation in a conventional 2D electron 
gas (2DEG) system, the plateau at zero conductivity is missing due to the appearance 
of the n = 0 LL at Ep = 0 |13511137] . While recent QHE experiments confirm the LL 
crossing, the plateau shift and the specific sequence of plateaus of the Hall conductivity 
in bilayer |13811140] and trilayer graphenes |14111146] . further explorations on the cases of 
distinct graphene systems are undertaken in order to comprehend the LL anticrossings and 
the unusually sequenced LLs. 

The magneto-electronic properties are evidently derived within the framework of the 
generalized tight-binding model, which is applicable to multilayer systems with various 
stacking conhgurations and layer numbers. The sublattice dominance is determined by the 
subenvelop functions spanned by the bases in the tight-binding model. Under a variation of 
the magnetic held strength, the LL spectra with respect to different stacking conhgurations 
can be realized. In terms of the relationship between the Landau state and the sublattice 
dominance, one can straightforwardly classify the LLs and dehne their quantum numbers. 
This provides a useful means for analyzing other essential properties of graphene systems. 
This review article broadly covers the helds related to multilayer graphene systems and 
discusses many factors ahecting their magneto-electronic properties. In Chapter 2, we in¬ 
troduce the Peierls tight-binding model to obtain the Hamiltonian matrices describing the 
AA, AB and ABC stacking conhgurations, in the presence of a perpendicular uniform mag¬ 
netic held. The Chapters 3-6 deal with the magneto-electronic properties of monolayer and 
multilayer graphenes with AA, AB and ABC stackings. We start with the case of mono- 
layer graphene in Chapter 3, and subsequently discuss AA-stacked graphenes in Chapter 
4, in which all layers are exactly aligned, leading to a symmetry-protected Dirac-like elec- 










tronic spectrum. The LL spectra of both monolayer and AA-stacked multilayer graphenes 
show the square-root energy dependence on the field strength and the quantum number. 
Chapters 5 and 6 are devoted to AB- and ABC-stacked graphenes, which both show a 
non-monotonic LL energy dependence, and reveal a complex pattern of LL anticrossings 
as a result of certain interlayer atomic interactions. In each of these chapters, we present 
the theoretical calculations of the LL spectra and the wave functions, certain important 
characteristics of which have been confirmed experimentally by STS, optical spectroscopy 
and quantum transport measurements (described in the final subchapter). We focus on 
magneto-absorbtion experiments, which are available for the investigation of the layer and 
stacking dependence of the magneto-electronic properties over a broad range of energies. 
Meanwhile, a detail comparison among these graphene systems is also presented in Chapter 
7. Finally, Chapter 8 contains concluding remarks. 


2 The generalized tight-binding model 

Electronic states are evolved to discrete dispersionless LLs in the presence of a magnetic 
held. In this chapter, we show how to calculate the Landau energy bands of multilayer 
graphenes in the framework of the generalized tight-binding model, which is based on the 
subenvelope functions of different sublattices. The Hamiltonian is built from the bases of 
tight-binding functions coupled with a Pierer phase factor. Three typical kinds of stacking 
conhgurations are used as study models: AA, AB and ABC. The magnetic quantization 
of electrons in graphene systems show interesting phenomena as a function the stacking 
conhguration and the number of graphene layers. 
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2.1 Monolayer graphene 


An illustration of the honeycomb structure of monolayer graphene is shown in Fig. 1(a). 
The primitive unit cell containing two carbon atoms, A (black) and B (red) atoms, is 
marked by the gray shadow. The symbols ao (— —2.6 eV) and }) (= 1.42 A) indicate the 
nearest-neighbor hopping integral and C-C bond length, respectively. Furthermore, ao is 
the most important atomic interaction often introduced to calculate the physical properties 
of monolayer graphene. The Brillouin zone corresponding to the primitive unit cell is shown 
in Fig. 1(b), where F, M and K are three highly symmetric points. Based on the primitive 
unit cell shown in Fig. 1(a), the wave function is expressed as T = where pa 

and pb, respectively, stand for the Bloch wave functions of the A and B atoms and are 
represented as USB] 


Pa = Aexp(*k ■ RA)x(r - Ra) 


( 1 ) 


Pb = Rexp(2k ■ RB)x(r - Rg). 

x(r — R^) (x(r — Rs)) is the normalized orbital 2pz wave function for an isolated atom at 
R^ (Rs), with A = ±1 indicating the bonding and anti-bonding forms. 

In the presence of a uniform perpendicular magnetic held B = Bq'z, the vector potential 
is chosen as A{r)=Boxy and the related periodic Peierls phaseAG^m' = ^ "^(r) • dr 

is introduced in the tight-binding functions, where 0o = hc/e (4.1356 x 10“^^ [T m^]) is 
the hux quantum [731475] . Owing to the periodicity of the Peierls phase, the primitive unit 
cell becomes an enlarged rectangle unit cell along the x-direction (armchair direction), as 
indicated in Fig. 1(c), where Rb is associated with the period along x and dehned as Rb = 
~ Accordingly, the rectangle cell includes ARb atoms {2Rb A and 2Rb 

B atoms); its length along the x-direction is given by I = This implies that the Bloch 
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wave functions under a uniform magnetic field can be expressed by the linear superposition 
of the ARb Peierls tight-binding functions in the rectangular unit cell: |Aik), l-Bik), |^ 2 k), 

|i? 2 k),.|A 2 ij-ik), |i? 2 /?-ik), |^ 2 ijk); \B 2 r\,)- To solve the ARb x ARb Hamiltonian matrix 

more efficiently, a band-like Hamiltonian matrix [nSlESHZS] is introduced by rearranging 
the bases as the sequence |Hik), \B 2 m), |Hik), \A 2 Rk), |H 2 k), |H 2 H-ik), |H 2 k), \A 2 R-ik), 

.I^ij-ik), \Br+ 2 \i), \BR-ik), |H/j+ 2 k), l^ijk), |H/{+ik), |H^k); I^^R+ik)- This matrix is 

thus expressed as 



0 

q* 

P*i 

0 


0 

0 

0 


Q 

0 

0 

P2R 

0 


0 

0 


Pi 

0 

0 

0 




0 


0 

P*2R 

0 

0 



0 




0 




0 

P*R 

0 


0 




0 

0 

0 

Pr+i 


0 

0 


0 

PR 

0 

0 

Q 


0 

0 

0 


0 

P*R+1 

q* 

0 

The elements pm and q are 

dehned as 







( 2 ) 


Pm = 7o exp -i [{ka:b/2 -F kyy/3b/2 + AGmm') + {kxb/2 - ky^/^b/2 - AGmm')] , 

q ='fo exp[i{kxb)]-, (3) 

AGmm' = -vr(m - 5/6)/R. 

By diagonalizing the Hamiltonian matrix, the energy dispersion and the wave function 
^c,!; obtained. 


11 






2.2 Multilayer graphenes 


When graphene layers are stacked along the 2 ;-direction, the physical properties of an N- 
layer graphene are strongly affected by its layer number, interlayer atomic interactions and 
stacking configurations [89ll^fT2lbl271fT3llll46] . Three typical stacking configurations, 
namely AA, AB and ABC, are usually chosen for theoretical and experimental studies. 
They present distinct physical properties, e.g., optical absorption spectra [HSHHI], Coulomb 
excitations |12411127] and the quantum Hall effects |13411146] . which do reflect the influences 
of layer number, stacking configurations and interlayer interactions. In this chapter, the 
Hamiltonian matrix of the AA-, AB- and ABC-stacked A^-layer graphenes are derived in 
the following. 

2.2.1 A A stacking 

The geometric structure of an AA-stacked A-layer graphene is illustrated on the top 
right side in Fig. 2. With an AA-stacking configuration, the geometric structure is 
formed by periodically stacked monolayer graphenes along the 5 direction. Four im¬ 
portant atomic interactions are included in the tight-binding calculations: one nearest- 
neighbor intralayer interaction ao = 2.569 eV, and three interlayer interactions ai = 
0.361eV, 02 = 0.013 eV and = —0.032 eV |209j . The primitive unit cell under a 
uniform magnetic field consists of A x ARb atoms. In order to obtain a band-like ma¬ 
trix, we arrange the bases as follows: \A\^), l^ik)> l^ik)l^rL), 

l-^lk); . I^2k)) |-®2k)) l^ik); |-®2k)vl^^ 1-^^ ^), |A^), . I^^sk); 

l^ 2 ^sk)> I^^Rsk); \B 2 Rsk)- One thing that should be noted 
is that all hopping terms and external fields considered in calculations regarding A-layer 
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graphenes are not just treated as perturbation terms; rather, the effects of atomic inter¬ 
actions and external fields can be simultaneously included in the calculations. For an 
AA-stacked TV-layer graphene, the nonzero matrix elements are 

I FT I Cyoitlk^m' ,m “1“ i4k^m' ^m',m “1“ ,m—l')^L—l,L' 

(^mkl-^l^m'k) = (^mkl-^l^m'k) = + a2(5L,L'+2, 

where the superscript in the tight-binding function represents the Lth layer. The four 

independent phase terms are 

tik = exp{i[—{kxb/2) — {^/3kyb/2)] + 7r(m — 5/6)/i?} 

+exp{i[—{kxb/2) + {y/3kyb/2)] — 7r(m — 5/6)/i?}, 

t2k = exp{i[-{kxb/2) - {y/3kyb/2)] + 7r(m - 3/6)/R} 

+exp{i[—{kxb/2) + {y/3kyb/2)] — 7r(m — 3/6)/i?}, (5) 

tsk = exp{i[-{kxb/2) - {y/3kyb/2)] 7r(m - l/6)/i?} 

-|-exp{i[— (fca;6/2) -|- {\/3kyb/2)] — 7r(m — l/6)/i?}, 

Uk = exp{ikxb). 

When an electric held is applied, the electric potentials are added to the site energies 
of the diagonal Hamiltonian matrix elements. 

2.2.2 AB stacking 

The diagrammatic drawing of an A-layer graphene in the AB-stacking conhguration is 
shown in Fig. 3. Each of the graphene sheets is obtained by shifting adjacent layers along 
the armchair direction by one bond length b' . The A atoms on any odd-number layer 
have the same (x, y) coordinates as B atoms on the upper layer, while its B atoms are 
projected at the center of the hexagon of the adjacent layer. The atomic interactions based 
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on the SWMcC model include one nearest-neighbor intralayer interaction 70 = —3.12 eV, 
five interlayer interactions, 71 = 0.38 eV, 72 = —0.021eV, 73 = 0.28 eV, 74 = 0.12 eV, 
75 = —0.003 eV, and the chemical environment difference between A and B jq = —0.0366 
eV IISIIIT]. The primitive unit cell of AB-stacked A^-layer graphene consists of 2 x 
atoms, and the magnetically enlarged rectangular cell comprises N x ARb atoms. Unlike 
the effective-mass calculation that treat non-vertical interactions as perturbation terms, 
such as 73 and 74 , the tight-binding calculations can well define the dimensions of the 
Hamiltonian matrix by simultaneously considering all the all SWMcC parameters and the 
magnetic field. The Hamiltonian matrix elements are expressed as 

(AjjjIhf I 'yi^m,mAL,L'±l T (T5 T L,L' T ~2^m,mAL,L'±2i 

— 73 (^ 2 ,+ tlkSm,m'-l)5L,L'±l + ^ L,L' + 

(AjyjIhf 15^/) 70(^4,fc^Tn,m' T + 1 ')^L,L' T 74(U,fc*^Tn,m' T fo^ odd T, 

7o(U,A:<5m,m' T +l)^L,L' T 74(^4,fc<^m,m' T +l)^L,L'±.1 fo'^ eveu L. 

2.2.3 ABC stacking 

The ABC stacking configuration is illustrated in Fig. 4, where each of the graphene sheets 
is shifted by a distance of b along the armchair direction with respect to the adjacent layer. 
The sublattice B of one layer is situated directly above the A atom of the adjacent lower 
layer, whereas sublattice A lies above the center of its hexagon. The calculation for such a 
system takes one nearest-neighbor intralayer interaction /3o = —3.16 eV and five interlayer 
interactions /9i = 0.36 eV, (52 = —0.01 eV, /33 = 0.32 eV, (3^ = 0.03 eV and (5^ = 0.0065 
eV, into account [16]. Consequently, the magnetically enlarged rectangular cell of an ABC- 
stacked A-layer graphene is composed of A x 4/?^ atoms and the Hamiltonian matrix, 
built by the ANRb bases, is represented as 


14 


{BLWA'i.y) 

— “t“ /^S(^3^^*5m,7n^+l “1“ L^V—1 ^ 3/ 2, 

(BLmA-i^y) 

/3o(^3,/i:*5m,m^+l H“ L^U H“ —l)*^L,L^—1 ^ 3/ 1, 

/3o(^2,A:*5m,m^ H“ ^4,/i:*5m,m' —l)^L,L^ H“ /33“1“ —l)^L,L^—1 3/, 

= /? 2 ( 5 m,m'< 5 L,L '-2 + I3i5m,m'SL,L '-1 for any L, 

{BL\B\BB^} = {Ai,f\H\AB*^) 


(^ 4 (^ 1 ,k^m,m‘ 

' + 

t4,k^m,m‘ 

1 

7 

for L = 

31- 

2, 


) = 

/ R^+11 
■ \-°mk 1 





/^4 

’+1 

+ t4,k^m 

,m')SL,L'-l 

for L = 

31- 

2, 

(Atk|77|Swk 

) = 

- 

- \-°mk 


— \^mk 

|B|4lti> 

/^5 (^3,/i:^m,7n' 

’+1 

T 

,m')^L,L'-2 

for L = 

31- 

2, 

(aLWab.^ 

) = 

= {-Blit 






' + 

tA,kSm,m‘ 

'-i)Sl,L'-1 

for L = 

31- 

1 , 

{BL\H\BB,^ 

) = 

-- -- 

- \^mk 

|B|Bli‘> 

/^5(j'2,kSm,m‘ 

’ + 

t4,k^m,m‘ 

'-i)Sl,L'+2 

for L = 

31, 


(Amk 

) = 

- 


— / A^'~^‘^ 

\^mk 



’ + 

t4,k^m,m‘ 

’-l)^L,L'+2 

for L = 

31; 



(7) 


The generalized tight-binding model can deal with the effects of several kinds of external 
fields in the calculations, e.g. modulated magnetic fields [83], periodic potentials [ 8 l| and 
composite fields [HS]. An exact diagonalization method that transforms the Hamiltonian 
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matrix into a band-like one is utilized to characterize the magnetic quantum numbers and 
wave functions. The effects of external fields and geometric structures are directly reflected 
in the electronic properties. The physical phenomena are explain in a clear physical picture. 
This model can be further developed to investigate magneto-optical properties and Coulomb 
excitations, which both are closely related to Landau-level spectra and transition matrix 
elements. The transition intensity needs to be evaluated from the initial- and hnal-state 
wave functions, in which each individual wave function is a linear combination of the 
products between the subenvelop function and the tight-binding function on each site in 
the magnetically enlarged unit cell. 

Generally, the generalized tight-binding model with an effective digitalization method 
is applicable to study the quantization effect in arbitrarily stacked graphenes, other layered 
materials, e.g. M 0 S 2 |114[|115ll222[|223j andsilicene j224fl226j . and bulk topological insulator 
|22711229] . The characteristics of electronic structures and wave functions are well depicted. 
Moreover, the results are accurate and reliable within a wide energy range. On the other 
hand, the effective-mass model expands the low-energy Hamiltonian near the K point, 
deriving an analogue of the Dirac Hamiltonian for massless fermions, where the interlayer 
atomic interactions and external helds are regarded as perturbations in the calculations. 
It could conceivably be used to comprehend the magnetic quantization at low energies for 
fewer-layer graphenes with symmetric structures. 
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3 Magnetic quantization in monolayer graphene 


3.1 Dirac-cone band structure 

Monolayer graphene exhibits the Dirac-type linear dispersion near Ep = 0 dne to the 
special hexagonal lattice strnctnre, as shown in Fig. 5. A magnetic held causes the states 
to congregate and induces dispersionless LLs, as shown in Fig. 6(a) for Bq = 20 T. The 
unoccupied LLs and occupied LLs are symmetric about Ep = 0, and characterized by 
quantum numbers which corresponds to the number of zeros in the wave function. 
Each Dirac-LL is fourfold degenerate without considering the spin degeneracy. Its energy 
can be described by the simple square-root relationship \E^'"\ oc 

3.2 Landau levels and Landau wave functions 

For understanding the magneto-electronic properties, researching the evolution of the LL 
wave functions is very useful. Based on the arrangement of odd-indexed and even-indexed 
atoms in the primitive cell (Fig. 1(c)), the wave function can be expressed as 

2Rb-1 2Rb 

|>i>k> = (ATM + BTM) + Y.(AT\M + Bp\b„^)), (8) 

m=l m=l 

where o(e) represents an odd (even) integer. Subenvelop functions of {A^) and 
{BE) stand for the probability amplitudes of the wave functions of the A and B atoms, 
respectively. The wave functions associated with the odd-indexed and even-indexed atoms 
in Eq. (7) have only a phase difference of vr, that is, AE = —AE and BE = —BE- 
Therefore, it is sufficient to comprehend the main characteristics of the LLs by discussing 
only the amplitudes of AE and BE- As shown in Figs. 6(b) and 6(c), the subenvelop 
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functions exhibit even and odd spatial symmetries and can be described by an order 

Hermite polynomial and Gaussian function. They are distributed around the localization 
center, that is, at the 5/6 position of the enlarged unit cell. On the other hand, similar 
localization centers corresponding to the other degenerate states occur at the 1/6, 2/6, and 
4/6 positions. A simple relationship of subenvelop functions exists between the A and B 
sublattices, i.e., of is linearly proportional to of + 1. Moreover, the con¬ 
duction and valence wave functions are related to each other by the following relationship: 
Aq = A^ and B‘^ = —B^. 

The dependences of the LL energies on and Bq behave differently at low and high 
energies, as shown in Figs. 7(a) and 7(b). The LLs with < O.dyo are proportional to 
the product of \/rF and to \/5o simultaneously, but there is no such relation for the LLs 
with > 0 . 470 , because the high-energy LLs are too densely packed to be separated from 
one another |64] . 

3.3 Oscillating Landau subbands under a modulated magnetic 
field 

As compared with the uniform case, a modulated magnetic held has a different impact 
on the electronic properties. The spatially periodic held along the armchair direction is 
assumed to be BMsm{27ix/3b'RM)'z, resulting in oscillating Landau subbands, as shown 
in Fig. 8 for the modulation held strength Bm = 20 T and the period Rm = 500. In 
the lower energy region, the conduction and valence subbands, symmetric about Ep = 0, 
appear around ky = 2/3. The subbands at Ep = 0 are partially hat. For the other Landau 
subbands with weak energy dispersions, they are characterized by double degeneracy and 
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have one band-edge state at ky = 2/3 and four extra band-edge states on both sides of 
ky = 2/3. The former band-edge state energies at ky = 2/3 are close to the LL energies, 
while the latter state energies are dependent on the modulation held strength and period. 

They demonstrate the strongest dispersion and the destruction of the double degeneracy. 
Moreover, the ky range with respect to both the weak dispersions and partial hat bands 
grows with an increasing held strength and a longer modulation period. On the contrary, 
when the inhuence of the modulation held becomes much weaker with increasing energy, 
the parabolic subbands in the higher energy region are similar to the twofold degenerate 
subbands directly obtained by the zone folding of monolayer graphene in the Bm = 0 case. 

3.4 Geometry-induced non-uniform magnetic field in curved graphene 
systems 

Investigations on curved systems (e.g., curved graphene nanoribbons and carbon nanotubes 
[22] - l2^[93] 1 provide another means of studying the ehect of a non-uniform magnetic held 
on graphene systems, and the cooperation and competition between the geometry and a 
magnetic held. In hnite-size ID graphene nanoribbons, the quantum conhnement limits the 
formation of Landau states in the low-energy region [261128113011 10 3] . where the energy bands 
are a combination of quaisi-LLs and parabolic bands, as shown in Fig. 9(a). The qausi- 
LLs possess the dispersionless energy parts in a constrained wave-vector range, whereas it 
increases in the cases of stronger magnetic hux, i.e., stronger helds or wider ribbons. As a 
hat ribbon is bent into a curved one, the quasi-LLs are depressed due to a weak magnetic 
quantization. As shown in Fig. 9(b), the Quasi-LLs entirely evolve into oscillating Landau 
subbands, of which the energy dispersion and the band-edge states are determined by the 
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confinements of both the open edge and the local magnetic held. On the other hand, 
in a carbon nanotube, a cylindrical structure corresponding to a zipped curved graphene 
ribbon, the zero net magnetic hux enhances the difficulty in forming the quasi-LLs. Instead, 
a transverse magnetic held results in a coupling of angular-momentum states, as shown in 
Fig. 9(c). With the magnetic held varying from low- to high- strength regimes, the coupling 
of independent angular momenta becomes signihcant. Recently in an extremely severe 
condition, i.e., the nanotube diameter is larger than the magnetic length, an experimental 
evidence of the Landau states in multi-walled carbon nanotubes under a very strong held, 
as high as 60 T, was identihed by magneto-transport measurements [SUES- 

Graphene, a hexible sheet of carbon atoms, can be curled or twisted without loosing 
its atomic structure and unique electronic properties. To date, a series of curled and 
folded graphene-related materials that posses specihc geometric symmetries and dimensions 
have been synthesized—such as graphene bubbles |210] , graphene ripples |211[|212j , carbon 
tori [in], carbon nanoscrolls ra, coiled carbon nanotubes |214H2l8] . Of further research 
interest are hrstly the quantization of electronic states resulting from the cooperation and 
competition between the geometry and a magnetic held, and secondly whether the quasi- 
Landau states can survive or not. 


4 Magnetic quantization in AA-stacked graphene 

AA-stacked graphene possesses a highly symmetric geometry structure, in which the carbon 
atoms of one layer are directly above/below the carbon atoms of another layer. In natural 
graphite, graphene layers are normally stacked in AB and ABC orders, while AA stacking 
conhguration is not found due to the much higher total energy in comparison with AB and 
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ABC configurations. Recently, however, it can be artificially synthesized and maintain the 


stability for macroscopically long times |1541I155] . The electronic band structure in the low- 
energy region comprises several Dirac cones [32]. The massless Dirac fermions with linear 
dispersion were verihed by STS experiments on the bilayer AA-stacked graphene [isa. 
The nature of the Dirac fermions attracts intense interest of theoretical and experimental 
research in AA-stacked graphenes. In this chapter, we mainly focus on how the low-lying 
electronic states within the Dirac cone are quantized by a uniform perpendicular magnetic 
held. Furthermore, we illustrate the evolution of the Landau-level spectra as a function of 
the total layers. The magneto-optical properties and the Coulomb excitations are closely 
related to the Landau-level spectra and the transition matrix elements that are evaluated 
from the initial- and hnal-state wave functions. 

4.1 Multiple Dirac-cone band structure 

In the absence of external helds, the A^-layer AA stacked graphene, most symmetric struc¬ 
ture among all the stackings, owns N sets of linear energy subbands intersecting at the 
K point, each behaving as the monolayer-like Dirac cone [32ll33l[66]. For the odd-layer 
scenario, one of the N sets is situated close to Ep = 0, mapped to the monolayer, and 
the other ones are oppositely mirrored to each other with respect to Ep = 0. For the 
even-layer scenario, half of the Dirac-cone structures are placed above Ep=0 and the an¬ 
other half below Ep = 0, the mirroring being in accordance with the odd-layer case. The 
trilayer case is illustrated to understand the main features of the energy dispersions. Three 
sets of subbands are shown in Fig. 10(a) and labeled by indexes S^’'^ and The 
Dirac cone structure of is nearest to Ep = 0, and those of and are approx- 
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imately symmetric to each other. Similarly, examples of even layers {N=2 and A^=4) 
are depicted in Fig. 10(b). It should be noted that the Dirac points at the K point are 
shifted relatively to their neighbor by the energy that decreases with the number of lay¬ 
ers, as well as demonstrated in the effective-mass model; the energy difference is given by 
2 cos(j -|- l 7 r/(A^ -|- l))-2 cos(j 7 r /{N + 1)) [ 66 ] . 

4.2 Landau levels and Landau wave functions 

As it is the general and fewest layer case, we select trilayer graphene to elaborate on the 
magnetic quantization of the Dirac quasiparticles in the AA stacking conhguration. The 
LLs at i?o = 20 T can be divided into three groups, as indicated by black, red and blue colors 
in Fig. 11(a), in which each Landau state at a given {kx,ky) is four-fold degenerate; the 
division is based on the energy spectrum and the characteristics of the spatial distribution 
of the wave function (Fig. 11(b)-11(g)). In particular, the conduction LLs and valence 
LLs of each group are approximately symmetric to each other about the zero-mode LL, 
= 0 , 77 . 2 ’^ = 0, or 77 , 3 ’' = 0 , as corresponding to the monolayer-like energy spectrum. 

The wave functions, characterized by the subenvelope functions on different sublattices, 
are used to precisely dehne the quantum number of each LL. For the four degenerate Landau 
states at {k^ = 0, = 0), the wave functions at the 1/6 and 4/6 locations of the enlarged 

cell have identical amplitudes, as do those at 2/6 and 5/6. Moreover, the 1/6 localized 
state has the subenvelope functions A(,’s, which are equal to iJ^’s (/ = 1, 2 and 3) of the 
2/6 localized state and vice versa. In Figs. 11(b)-11(g), the subenvelop functions of the 
2/6 localized state are demonstrated, where further discussions are restricted to the odd- 
indexed components because of the relationship With the starting energy away 
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from Ep = 0, the first and third groups of LLs have a monolayer-like relationship between 
the two subenvelop functions for every layer, while the second group of LLs spreading 
from Ep = 0 have the same relationship only in the two outmost layers. The subenvelop 
functions related to the inner layer are empty. This means that the LLs of the second group 
are entirely regarded as a result of monolayer graphenes, and the Landau quantization of 
the Dirac cones mapped to the monolayer is consistently presented in the odd-layer cases. 
On the contrary, all the layers of the even-layer AA-stacked graphenes are equivalent. The 
quantum numbers of LLs can be identihed by the non-empty subenvelop functions of lattice 
B in any layer. Using the single-mode wave functions is a visible means to identify the LLs, 
since there is an absence of LL couplings under the specihc interlayer atomic interactions. 
These results indeed reflect the monolayer properties of AA stacking. 

The unique spectral features are better understood by investigating the LL evolution 
with the held strength Bq. The three monolayer-like LL spectra of AA-stacked trilayer 
graphene are shown in Fig. 12. The energy distribution in each group LL is related to 
the individual Dirac cones that are primarily separated by the specihc interlayer atomic 
interaction ai [32l|33l|66]. Their energies can be regarded as those of the monolayer LLs, 
which have a simple square-root dependence on the held strength and the quantum number. 
Meanwhile, all the LLs directly cross one another, i.e., there is an absence of LL anticross¬ 
ings, since each Landau state is characterized by a single-mode harmonic function. 

The LL spectra of the AA-stacked bilayer and tetralayer graphenes are shown in Fig. 
13 and Fig. 14, which respectively demonstrate the two and four groups of monolayer-like 
patterns, as well as the oscillating Fermi level. Each group is roughly symmetric about the 
zero-mode LL, which is nearly constant as a function of held strength; each level moves 
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towards the individual Dirac point when the magnetic held approaches zero. There are 
essential differences between the tight-binding model and the effective-mass model. In 
systems with charge neutrality, the tight-binding model calculates the oscillating Fermi 
levels based on equivalent numbers of conduction and valence carriers. On the other hand, 
due to the Dirac cone structures extending to the full k space in the effective-mass model, 
the isotropic dispersion makes the spectra identical and the Fermi levels consistent with its 
zero-field value. 

4.3 Magnetic-field-dependent oscillation of the Fermi level 

Around the zero-field value of Ep = 0, the Fermi energy at a constant carrier density 
oscillates as a function of the field strength, as shown in Fig. 12 by the bold wriggling 
curves. This leads to an energy gap, the spacing between the highest occupied LL and the 
lowest unoccupied LL; this gap oscillates with the magnetic field, and results in the metal- 
semiconductor transition at certain specific fields. The magneto-resistance measurements 
to detect the oscillation of the Fermi level were performed in other 2D electron gas systems 
|21911221] . As for the i?o-dependent energy gaps in bilayer and trilayer graphenes (black 
curves in Figs. 15(a) and 15(b)), the oscillation period and amplitude both increase with 
an increment of the field strength, but decrease in the cases of more layers. Furthermore, 
a gate voltage shifts the Dirac cones [33], providing an available means of using external 
fields to modulate the energy gap. 
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4.4 Magneto absorption spectrum 


The main characteristics of the electronic properties is directly reflected in other physical 
properties, e.g., optical absorption spectra. In the absence of external helds, previous 
works show that only intra-subband excitations within the same Dirac-cone structure are 
permitted, and consequently asymmetric absorption peaks are induced [86]. The threshold 
excitation energies are investigated in the odd-layer and even-layer scenarios. Also shown 
are their dependence on the stacking layers and, after exceeding certain layer numbers, 
the cases in both scenarios exhibit critical properties identical to those of 3D graphite 
for different frequency regions. On the other hand, feature-rich spectra are expected due 
to the optical transitions between LLs. There are several kinds of prominent symmetric 
peaks, resulting from the individual intragroup LL transitions and satisfying the particular 
selection rule An = ±1 (same as the monolayer) [65]. Transitions between intergroup LLs 
are forbidden due to the zero electric dipole moments derived from the spatial symmetries 
of the wave functions in the AA-stacked systems. The effective-mass model demonstrates 
similar spectra [66] . However, as a result of the identical spectra of groups of LLs and the 
constant Fermi level, important characteristics were missed in the case of the AA stacking, 
such as the strength and numbers of peaks and the excitation channels. 

5 Magnetic quantization in AB-stacked graphene 

The AB stacking conhguration is stable and commonly observed in nature graphite. Gen¬ 
erally, the AB-stacked graphenes have attracted intense attention due to the fact that 
the properties of massless Dirac fermions are presented in all cases of an odd number 
of graphene layers [35h[39] . Of special interest is the band structure that can be mod- 
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eled by a hybridization of a monolayer and bilayer graphenes. On the other hand, in all 
cases of an even number of layers, the energy dispersion is parabolic in the low-energy 
region. An evidence of Dirac nature of charge carriers in ultrathin epitaxial graphenes 
was observed by scanning tunneling spectroscopy |162j and angle-resolved photoemission 
spectroscopy |1821I1841I185] . The Landau quantization gives rise to interesting phenomena 
of the electronic properties in the even-layer cases of inversion-symmetry graphenes and in 
the odd-layer cases of symmetry-broken graphenes. 

5.1 Band structure: hybridization with monolayer graphene 

The low-energy electronic properties in AB-stacked bilayer graphene are characterized by 
four parabolic energy bands, as shown in Fig 16(a): a pair of conduction and valence 
bands, labeled (black curves), and another, labeled S^’'^ (red curves). The hrst pair 
of subbands slightly overlaps near Ep = 0 (inset of Fig. 16(a)), and each subband of 
the second pair is approximately shifted away from Ep = 0 by an energy of 71 or — 71 , 
where 71 is the vertical nearest interlayer atomic interaction. Unlike only parabolic bands 
in the bilayer graphene, the band structure in a trilayer graphene is hybridized from a 
monolayer and bilayer band structures, as shown in Fig 16(b). Near the K point, the 
intersections of linear and parabolic subbands are, respectively, separated, but a slight 
overlap remains between the valence and conduction bands (inset of Fig. 16(b)). All 
odd-layer scenarios still exhibit monolayer-like linear dispersions, while all even-layer ones 
only show parabolic energy dispersions [3SH3H]- AB-stacked graphenes are gapless 2 D 
semimetals, whereas a tunable energy gap can be induced by applying an electric held in 
the bilayer case |45] - [55] . The theoretical results are experimentally conhrmed by infrared 
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spectroscopy [SSlEn], ARPES isa and transport |5B] - I5U] measurements. 

5.2 Landau levels and Landau wave functions 

The trilayer graphene is used to elaborate on the main characteristics of the magnetic 
quantization of AB-stacked graphenes. The LLs at Rq = 20 T can be divided into three 
groups, as indicated by black, red and blue colors in Fig. 17(a). The first, second and third 
groups are, respectively, attributed to the quantization of the linear subband Si^ and the 
parabolic subbands S^’'^ and The LL spacing is smaller in the higher group than in 
the lower group, a spacing variation determined by the zero-held energy dispersion (DOS). 
In the hrst group, the onset energy is close to Ep = 0, and each Landau state has the wave 
function distributed only on the two outer layers (black curves), which is attributed to the 
magnetic quantization effect on the linear subbands As a result of the degeneracy 

breaking of the Dirac points (inset of Fig. 16(b)), two double degenerate states are split 
from the zero mode of the quantized Dirac quasiparticles, with the respective quantum 
numbers = 0 (5 meV) and = 0 (—5 meV), being determined by the dominant 
subenvelope functions of (localized at 1/6 and 4/6 unit cell) and (localized at 
2/6 and 5/6 unit cell), as shown by the black curves in Figs. 17(b)-17(m). The higher 
conduction LLs are placed with narrower spacings. Their quantum numbers n1 are in an 
ascending order with an increase of energy, and can be derived from either or R/’^ as 
well. 

The second group has two sets of double degenerate LLs, the dominating sublattices of 
which are R^ for the 1/6 (4/6) localized states and R/’^ for the 2/6 (5/6) localized states, 
as shown by the red curves in Figs. 17(b)-17(m). The first four double degenerate states. 
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labeled n 2 = 0 , 1 and ^ 2 = 0 , 1, are split from the zero and first modes of the S 2 '’ subband 
quantization, with the splitting energies ~ 10 — 20 meV. However, the splitting energies 
for higher LLs are less than 1 meV, which might be difficult to observe in experimental 
measurements. |141H144] In the same way, the dominating sublattices are for the 1/6 & 
4/6 localized states, and for the 2/6 & 5/6 localized states (blue curves in Figs. 17(b)- 
17(m)) for the third group. In all, each sublattice equally dominates the magneto-electronic 
properties. 

5.3 Symmetry-breaking effect on Landau level spectrum 

The LL spectrum of the trilayer graphene demonstrates the monolayer-like LLs in the hrst 
group and the bilayer-like LLs in the second and third groups [691172] . as shown in Fig. 
18. The onset energies of the conduction LLs in the hrst, second and third groups are 
identihed as the band-edge state energies of subbands SI, S '2 and S'g, respectively, and the 
same identihcation is applied to the valence LLs. For the hrst group, the energy diherence 
between the onset valence and conduction LLs is equal to the energy splitting of the Dirac 
points at the K point, and for the second group it is determined by the gap between the 
S '2 and S'! parabolic subbands. Being split from Ep = {], the 12 (4iV) Dirac-point LLs, 
including the two double degenerate ones of n^’ = 0 and the four double degenerate ones 
of — 0 conhned in a narrow energy range (~ 10 meV) and insensitive 

to the held strength. Besides, the splittings between the 1/6 (4/6) and 2/6 (5/6) localized 
states are revealed for the low-lying LLs. The intergroup LLs frequently cross one another, 
but the anticrossings between them only occasionally appear in strong helds. The quantum 
transport experiments have verihed the degeneracy change of the low-lying LLs due to the 




lack of inversion symmetry and the LL crossings in AB-stacked trilayer graphenes m- 

Essentially, the LL spectra of AB-stacked graphenes can be classihed into two kinds 
based on whether the numbers of their stacking layers are even or odd [67f[72] . To com¬ 
prehend the general results, we illustrate the spectra of even-layer AB-stacked graphenes. 
Unlike the degeneracy breaking in the trilayer graphene, the LL spectra of even-layer AB- 
stacked graphenes all consist of four-fold degenerate LLs due to the inversion symmetry of 
the geometric structures, such as the spectra of bilayer and tetralayer graphenes shown in 
Fig. 19 and Fig. 20, respectively. The held-dependent energies depend on the subband 
dispersions in a zero held. One can End the dominant sublattice for each group, i.e. in the 
ascending sequence of groups followed by (B^) and Al (A^) for bilayer graphene [ 68 ] . 
and Bl (B^), B^ (B^), Al (A^), and A^ (A^) for tetralayer graphene [75]. The 4A^ Dirac- 
point related LLs, insensitive to magnetic helds, are conhned in the vicinity of ^^^=0, and 
classihed as n'l=0 and nl= 0 for bilayer graphene, or 71^=0, n\=l, '^■ 2 = 0 , and ^ 2=1 for 
tetralayer graphene. 

As a result, it is deduced that for the LLs in A^-layer graphene, there are [(A^ -|- l)/2] 
groups starting to appear near Ep = 0, of which one especially exhibits a monolayer-like 
energy spectrum once N is odd, while the other N — [(A^ -|- l)/2] groups start to appear 
away from Ep = 0. The degeneracy of the LLs is broken in the odd-layer cases, while 
it is preserved in the even-layer cases due to the inversion symmetry. Furthermore, the 
anticrossing of LLs between two neighboring groups occurs in the tetralayer graphene at 
certain regions of magnetic held strengths and energies, as the ellipses depict in Fig. 20. It 
is revealed that the interlayer atomic interaction 73 couples the LLs in the 3 modulo [HSl 
[7811208] , so that according to the Wigner-von Neumann non-crossing rule, two multi-mode 
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LLs containing identical modes do not cross each other as a function of the magnetic held 
strength. An explanation is made in detail in the following chapter 6.3. With an increasing 
number of layers, the LL anticrossings are easily observed. While experimental observations 
on QHE conhrm the low-energy spectra through the specihc sequence of Hall conductivity 
plateaus in the AB-stacked bilayer [13911140] and trilayer graphenes [1411114511146] . their 
dependences in multilayer systems need to be systematically verihed. 

The dominating sublattices for each subband under magnetic quantization obey a reg¬ 
ular rule. For odd-layer graphenes, the dominating sublattices for the hrst group of LLs 
quantized from the linear subbands are the two outmost layer sublattices A^, or 

. The other (A^ — groups starting from (away from) Ep = 0 are followed 

by sublattices B (A) being counted from the inner (outer) to the outer (inner) layers in the 
ascending sequence of groups. In the same way, even-layer graphenes own a rule similar 
to the former, but excluding the linear-band part. This generalization for the magneto- 
electronic properties is valuable in understanding other physical properties, such as the 
mechanisms of magneto-optical properties and Coulomb excitations, etc. 

5.4 Magneto absorption spectrum 

The optical absorption spectra reflect the main characteristics of the electronic properties. 
In the absence of helds, only inter-subband excitations between the same pairs of parabolic 
subbands induce asymmetric absorption peaks [87] ; those excitations can be enhanced by 
applying an electric held [92] • A magnetic held, on the other hand, induces symmetric 
prominent delta-function-like peaks as a result of the intra- and inter-group LL transitions 
[97] . According to the Fermi’s golden rule, the former transitions satishes the monolayer 
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selection rule An = ±1, while the latter follows more complicated selection rules, which 
are not unihed for all inter-group transitions and should be identihed by checking the 
probability of the dipole transition. In a previous work |97] , an AB-stacked bilayer system is 
shown to exhibit an absorption spectrum consisting of four types of absorption peaks, with 
two being the result of intra-group transitions and the other two being the result from inter¬ 
group transitions. Each type of optical excitations lead to twin-peak structure owing to the 
asymmetry of valence and conduction LLs about the Fermi level (electron-hole asymmetry). 
The inter-group selection rules are Ani 2 = 0, Ani 2 = —2, AU 21 = 0, and AU 21 = 2 where 
the subscript index represents the group number. It is evident that the two different rules 
are both based on the underlying concept that the initial and hnal states in the dipole 
transition have the same zero points in subenvelop functions. This concept provides an 
insight into multilayer systems. The specihc selection rules for any inter-group transitions 
can be straightforwardly derived from the wave functions. It should be noted that for the 
odd-layer cases, the excitations from the quantized LLs of the linear bands to those of 
the parabolic bands are forbidden due to the zero dipole transition moment. Therefore, 
A-layer AB-stacked graphene is deduced to have a total of types of absorption peaks, 
resulting from A intra- and (A^ — A) inter-group transitions for even-layer cases; a total 
of ((A — 1)^ -f 1) types including A intra- and ((A — 1)^ — (A — 1)) inter-group transitions 
for odd-layer cases. The twin-peak structures are presented in recent magneto-absorption 
experiments [TIMfTTn] . where the measured transition energies, however, are limited to few 
hundreds of meV. The magneto-optical properties need to be verihed in detail over the 
entire spectrum range within the consideration of intra- and inter-group LL transitions. 
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5.5 Magneto plasmons 


In addition to optical properties, electronic excitations are also enriched by the magnetic 
quantization. The complicated Coulomb screening in multilayer graphene can be charac¬ 
terized by the dielectric-function matrix |1221I123] . The well-developed layer-dependent 
matrix is based on the tight-binding functions of distinct sublattices, as done for the mag¬ 
netic Hamiltonian matrix in the generalized tight-binding model. As a result, the interlayer 
e-e interactions, the interlayer atomic interactions and the magnetic held are simultaneously 
included in the calculations of Coulomb excitation spectra. By utilizing this matrix, the 
screened energy-loss function is obtained, which is then used to determine the main features 
of collective Coulomb excitations. The previous study shows that monolayer graphene ex¬ 
hibits a lot of low-frequency magneto-plasmons at zero temperature even if it is a zero-gap 
semiconductor nzg. Each magneto-plasmon mode is associated with a specihc inter-LL 
excitation channel from an occupied LL to another unoccupied one. Magneto-plasmons, 
quanta of the electron density oscillations in a perpendicular magnetic held, present the 
novel dependence in the relationship between their frequencies and the transferred momen¬ 
tum, since the longitudinal Coulomb force strongly competes with the transverse Lorentz 
force. The intensity, frequency, and number of magneto-plasmons are signihcantly ahected 
by the transferred momentum, the magnetic-held strength, the temperature and the doping 
carrier density |128fll33j . 

Such characteristics are further diversihed by the distinct stacking conhgurations. For 
AB-stacked and AA-stacked bilayer graphenes, there exist complicated energy-loss excita¬ 
tion spectra with more inter-LL plasmon resonance peaks |127j . In the former the energy- 
loss spectrum is mainly dominated by discrete inter-LL excitations, while in the latter a 
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2D-like plasmon involving the entire low-frequency Landau states is found. Also noticed 
that a 2D-like plasmon, being similar to that of a 2D electron gas in the absence of mag¬ 
netic field, is a result of the highly symmetric AA stacking and the dense LL distribution 
around the Fermi level. The inter-LL plasmon and the 2D-hke plasmon have a very differ¬ 
ent dependence on the transferred momentum and the magnetic held strength. These two 
kinds of magneto-plasmons are expected to be drastically changed by the increasing layer 
number. Whether there are new collective excitation modes in ABC-stacked graphenes 
deserves a further detailed study. 

6 Magnetic quantization in ABC-stacked graphene 

The ABC stacking conhguration stably exists in nature graphite jl891ll9U] . The surface- 
localized states of ABC-stacked graphenes form partially hat bands near the Fermi level 
[^3] - H5] . giving rise to an enhancement of the density of states of LLs. A prominent peak 
appears at the Fermi level with the intensity proportional to the layer numbers as a result 
of the increasing LLs that originate from the Dirac points. In particular, the band struc¬ 
ture exhibits sombrero-shaped energy dispersion around the energies of the nearest vertical 
intralayer hopping integral [391 - I12] . Consequently, the Landau-level spectrum reveals in¬ 
triguing features within such a band width ca: LLs with energies inversely proportional 
to the magnetic held strength and intra-group Landau-level anticrossings in the region of 
weak magnetic helds. 
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6.1 Band structure: surface-localized states and sombrero-shaped 


energy bands 

The low-energy electronic properties of ABC-stacked trilayer graphene are characterized 
by three pairs of conduction and valence subbands, S 2 '’ and Agdisplaying partially 
flat, sombrero-shaped and linear shaped dispersions, respectively, as shown in Fig. 21(a). 
In the vicinity of Ep = 0, the subbands of the hrst pair are weakly dispersive near the 
K point. As a result of the trigonal warping effect [73l[76], the low-energy dispersions are 
anisotropic with three crossing points at Ep = 0, which are identihed as Dirac points [75] . 
i.e., one in the KT direction and another two rotated with respect to the former by 120° 
and —120° (inset of Fig. 21(a)). The states in this region are formed predominantly by 
atomic orbitals on the two outermost layers, known as surface localized states [151 - H5] . 
Recent experiments confirm that a tunable energy gap can be induced by applying an 
electric held in the trilayer cases [5ll[58ll60]. Away from Ep = 0, the conduction subbands 
and the valence subbands, respectively, cross near the energies of (3i (=0.29 eV) and —f^i 
(=-0.29 eV), corresponding to the band edges of the subbands 82 ^ and In particular, 
the sombrero-shaped conduction (valence) subband, Af { 82 ), has a local energy minimum 
(maximum) and maximum (minimum), measured respectively, as 0.29 (-0.29) eV and 0.36 (- 
0.36) eV, within which two closed circular loops at a given energy lead to different quantized 
Landau orbitals. 

Tetralayer stacked graphene has a similar band structure, except for one more pair of 
parabolic bands located near ±/3i, as shown in Fig. 21(b). Surface localized states are 
distributed around Ep = 0 while four Dirac points come into existence, as shown in the 
inset of Fig. 21(b). The theoretical band structure is experimentally conhrmed by infrared 
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absorption spectroscopy jlU 2 j . For all ABC-stacked graphenes, the energy dispersion of 
the surface localized states is simply described by a relationship E oi. in the minimal 
model |39lll2lll3l[76] . Furthermore, it is shown that the number and characteristics of the 
Dirac points are related to the number of the stacking layers [SnilTSlES] • In addition, more 
of the subbands intersecting near ±71 show a sombrero-shaped energy dispersion with the 
increasing number of layers. The feature-rich band structure leads to the intricate and 
versatile magnetic quantization in ABC-stacked graphenes. 

6.2 Landau levels and Landau wave functions 

The magnetically quantized LLs of ABC-stacked trilayer graphene subjected to i?o = 25 T 
are shown in Fig. 22(a). The LLs can be divided into three groups according to the spatial 
characteristics of the wave functions, as indicated by black, red, and blue colors in Figs. 
22(b)-22(g). The hrst, second and third groups, respectively, result from the quantization 
of the weakly dispersive , the sombrero-shaped and the linear S' 3 ’’^ subbands. In 
particular, the onset energies and level spacings for all groups are relevant to the zero- 
held energy dispersions. Each level is four-fold degenerate due to the inversion symmetry, 
leading to the 1/6 (4/6) localized state having the subenvelop functions A(,’s (/ = 1,2 and 
3 for layer numbers) equal to —i?o’s (/ = 3, 2 and 1) of the 2/6 (5/6) localized state and 
vice versa. 

The hrst group starts to form near Ep = 0, in the vicinity of which the three degenerate 
LLs in the minimal model are slightly split by remote interlayer hoppings [731l75ll76j . Away 
from Ep = 0 the spacing becomes more obvious and the LLs are clustered together with 
increasing energy due to the higher DOS of the subband. As evidenced by the dominating 
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subenvelop function 5^ in Fig. 22 (c), the LLs are counted by increasing quantum numbers 
with increasing energies. On the contrary, the quantization related to the opposite outmost 
sites of is introduced to the 1/6 and 4/6 localized states. Near Pi, the LLs are densely 
packed as a consequence of the crossover of the and S'g subbands, and the second 
and third group LLs appear in the sequence of the quantum numbers determined by the 
subenvelop functions and in Figs. 22 (g) and 22(e). 

6.3 Intra- and inter-group Landau-level anticrossings in the field- 
dependent energy spectrum 

The three groups of LLs display feature-rich energy spectra during a variation of held 
strength, as shown in Fig. 23. Each LL moves towards the respective onset energy for 
every group, i.e., the local extreme value near the K point (Fig. 21(a)) with respect to each 
subband in the absence of an external held. The three Dirac-point related LLs are insensi¬ 
tive to the variation of held strength and conhned in a small energy range near the Fermi 
level, with their separated energies mainly determined by the specihc interlayer atomic 
interaction 74 . In contrast to what the monotonically B-dependent energies present in the 
hrst and third groups, a complicated pattern of non-monotonic LL dispersions is observed 
in the second group under certain critical helds. This is responsible for the particular 
magnetic quantization rule relating to the sombrero-shaped energy dispersion. There are 
several inverted LLs located between energies of the local maximum and minimum of this 
sombrero-shaped subband. The relationship between each LL and the critical magnetic 
held is determined by the interplay between the magnetic degeneracy and the total number 
of states enclosed in this area. A more detailed discussion is made as follows. 
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The LL spectrum exhibits special phenomena, the intragroup and intergroup LL an¬ 


ticrossings, the exploration of which is performed by following the evolution of the wave 
functions. It is seen that two LLs repel each other, as in the region of the dashed ellipses 
in Fig. 24, when these two levels are in modulo 3 [5^[75| 15211208] . This event can even be 
triggered among all different groups in a sufficiently strong magnetic held. This is a result 
of the Wigner-von Neumann non-crossing rule, which states that two multi-mode LLs do 
not cross as a function of the magnetic held strength if they simultaneously have identical 
modes. From a qualitative perturbation analysis, it is deduced that there are three factors, 
{32, ^3 and /^s, that induce the LL anticrossings in ABC-stacked graphenes. The energy 
perturbation caused by ^2 has its expectation value between the two same mode subenvelop 
functions and of the two hybridized states. A likewise explanation of the /Sa-induced 
correction is applicable to (3^, which is also responsible for the specihc relationship between 
and B^, and between and A^. The drastic changes of the wave functions during the 
anticrossing processes, e.g., for = 2 and Ug = 0 LLs, are illustrated in Figs. 24(b)-24(g). 
In principle, for Landau states away from the anticrossing regions (the hrst hve circles Fig. 
24(a)), the two subenvelop functions established on the perpendicular projected sublattices 
in adjacent layers have the same number of the zero points; two sub lattices in the same 
layer have zero-point numbers differed by one. As shown in the lower part of Fig. 24(b)- 
24(g), the non-perturbed = 2 LL has well-dehned subenvelop functions that have zero 
points of 3, 4, 2, 3, 1 and 2 in the order of {Al, B^, A^, Bl, Al, B^}. However, at the center 
of the anticrossing region, depicted in the middle part of Figs. 24(b)-24(g) (~0.37-0.41 
eV), the perturbed subenvelop functions possess multi modes composed of the main mode 
77-2 = 2 and the side mode rig = 0. This signihcant hybridization results in the two com- 
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ponents being comparable, while away from this region, the wave function converts to the 
characteristics of the well-defined rig = 0 Landau state (upper part of Figs. 24(b)-24(g)) 
as a result of the decreasing perturbation. 

A further exploration of the different layer cases is worthwhile in order to comprehend 
the unusually sequenced LLs and the LL anticrossings. Except for one pair of bands 
contributed by the surface localized states that is always weakly dispersive near the Fermi 
level, there are — 1 pairs of conduction and valence subbands intersecting in the vicinity 
of the energies ±/3i [H]. Some of them display the sombrero-shaped energy dispersions 
which determine the unique pattern in the LL spectrum. Tetralayer graphene shows only 
one group of LLs inversely dependent on the held strength in Fig. 25, and pentalayer 
graphene shows two such groups in Fig. 26. 

However, a relatively narrow region of the inverted LLs is observed for the graphene 
with fewer layers, which means that the enclosed states in the sombrero-shaped band 
structure are relatively small, as are the corresponding critical magnetic helds Bg. It is 
deduced that with an increase of stacking layers, Bg becomes stronger and the intragroup LL 
anticrossings show up in a relatively large energy range as a result of more sombrero-shaped 
subbands m- The derived rules are that the number of groups showing self anticrossing is, 
respectively, {N — l)/2 and {N — 2)/2 for odd and even N. Moreover, the generalization of 
the Dirac-point related LLs separated near Ep = 0 is also extended to the iV-layer systems, 
as shown in Fig. 27. There exist N four-fold degenerate LLs with the energy distribution 
insensitive to the magnetic field and limited in a very narrow region evaluated to be ~ 10 
meV for N > 10. The aforementioned results reveal feature-rich spectra in the cases of 
more layers. 


38 


6.4 Peculiar quantization rule for sombrero-shaped energy bands 


6.4.1 Inverse dependence of Landau-level energies on the magnetic field 

As a signature of the ABC-stacked graphenes, the intragroup LL anticrossings frequently 
appear in an energy region within the local maximum and the local minimum of the 
sombrero-shaped subband, as shown in Fig. 28(a), especially for weak fields where the 
inverted LLs is densely sequenced. We consider the Hamiltonian without the anticrossing 
interaction in order to clearly realize the particular quantization rule for the sombrero¬ 
shaped subband. The inverted curvature of such a subband structure gives rise to two 
available Landau orbitals at a given energy, one located at the inner circular loop and the 
other at the outer circular loop of the sombrero-shaped subband surface. A weak magnetic 
field aggregates only the states near the vertex of the inner surface, so that the LL energies 
at first decrease with increasing field strength as more inner-part states are quantized into 
LLs. At a critical magnetic field, labeled Bs{ns), the LL with the quantum number 
reaches the minimum of its energy distributions, which is related to the energy minimum 
of the sombrero-shaped subband, and then converts back to the typical i?o-dependence 
with further increasing field strength. This means that under the critical field Bg{ns), all 
the electronic states covered by the inner sombrero-shaped surface are magnetically quan¬ 
tized to the inverted LLs, in which the LLs are entirely filled up to the quantum number 
Ug- In contrast, the quantization of the states on the outer surface is introduced to the 
i?o-proportional LLs with quantum numbers larger than Ug. 
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6.4.2 Periodical quantization effect related to the inverse of the critical mag¬ 
netic field 

Given that the magnetic degeneracy of each LL is D = (eL^/27rhc) = pB, one can derive 
the relation Ng = Ug x D ai the critical point, where Ng is the nnmber of total states on the 
inner snrface. At a stronger magnetic held, fewer LLs can occnpy the totally enclosed states 
dne to the higher degeneracy D, and therefore the nnmber of the nnnsnal LLs decreases. 
Once the magnetic held is stronger than Bgirig = 1), a series of all the LLs evolves into an 
ascending seqnence with increasing energy. The relationship between Bg and Ug is simply 
characterized by Bg x Ug = const, as indicated by the constant A (l/i?o) period in Fig. 
28(b). Therefore the total degeneracy D x Ug eqnals the states enclosed in this area, i.e., 
Ng = D X Ug (X Bg X Ug, that is to say, the inversely proportional relationship between Bg 
and rig characterizes the magneto-electronic properties resulting the sombrero-shaped band 
structure. 

6.5 Magneto absorption spectrum 

In the absence of external helds, optical absorption spectra exhibit shoulder and peak 
structures in the low-frequency spectrum, dominated by the excitations from the surface 
localized states [SnilHI]- The absorption frequencies and intensities can be modulated by 
an applied electric held [50]. On the other hand, symmetric delta-function-like peaks are 
expected to be presented in the magneto absorption spectra as a result of the direct intra- 
and inter-group LL transitions. The particular selection rule An = ±1 is applied to the 
individual intragroup LL transitions. However, alike to the manner of the AB stacking, the 
allowed inter-group LL transitions are based on the underlying concept that in the dipole 
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transition, the initial and final states have the same zero points in the snbenvelop functions. 

For an iV-layer system, there are types of absorption peaks deduced from the transitions 
of which a total of N are from intra groups and a total of (iV^ — N) are from inter groups. 

Another interesting aspect of the absorption spectrum is that the LL anticrossing effect 
induces extra selection rules of modulo 3 for the intra-group and inter-group transitions. 

This is especially important for the intra-second group transitions under magnetic fields 
smaller than the critical strengths ca, where the frequently appearing anticrossings lead 
to a considerable number of peaks with intensities proportional to the LL hybridization. 

7 Differences among AA-, AB- and ABC-stacked graphenes 

7.1 Landau-level spectrum 

The stacking configuration plays an important role in the electronic properties of multilayer 
graphenes. Of interest is the respective energy band structures and the corresponding LL 
spectra of different stacking configurations, i.e., AA, AB and ABC stackings, as reviewed in 
the Chapters. 3-6. In the absence of external fields, the AA-stacked graphenes preserve the 
conical Dirac cone structures due to the vertically projected geometric structure. On the 
contrary, the AB- and ABC-stacked graphenes display rather distinct band structures, with 
significant hybridization caused by the interlayer atomic interactions. The band structure 
of AB-stacked graphene consists of several bilayer-like parabolic bands, while particularly 
in odd-layer cases, it owns additional monolayer-like linear bands. In ABC systems, there 
is always one pair of weakly dispersive energy bands near the Fermi level; those bands are 
formed by surface localized states. In addition, some bands peculiarly display sombrero- 
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shaped dispersions near the energy of the vertical nearest-neighbor interlayer interaction. 
The number of the sombrero-shaped subbands is increased when the number of layers is 
increased. For both AB- and ABC-stacked graphenes, several Dirac points are located at 
low energies with their characteristics related to the stacking conhguration and the layer 
number. 

The zero-held band structures determine the main characteristics of the quantized LL 
states, such as the sublattice dominance in Landau states, the held-dependent LL energies 
and the inter- and intra-group LL anticrossings. Each group of LLs is responsible for the 
magnetic quantization of the respective subband. The i?o-dependent energies and the onset 
energies of the LLs are determined by the energy dispersions and the local extreme values of 
the subband. According to the sublattice dominance for each group, the quantum numbers 
are well dehned based on the traits of the Landau wave functions. The LLs are four-fold 
degenerate in both AA and ABC systems, while the degeneracy is broken in odd-layer AB 
systems due to the lack of inversion symmetry. The broken degeneracy is demonstrated 
by the Landau wave functions with different localization centers. The i?o-insensitive LLs 
are presented for all the three stacking conhgurations, reflecting the properties of the Dirac 
points. However, unlike the large energy separations among them in AA systems, the 
distribution of such LLs is smaller than 10 meV for AB and ABC systems, with the slight 
LL splittings mainly caused by the interlayer atomic interaction [3^ ( 74 ). 

The intra- and inter-group LL anticrossings are important in detailing the magneto- 
electronic properties of multilayer graphenes. This can be attributed to the Wigner-von 
Neumann non-crossing rule which prevents the coupled LLs in the modulo 3 from crossing 
one another during a variation of the held strength. The AA-stacked graphenes present 
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direct crossings by either the same or different groups since there exist no couplings among 
LLs so that each LL is described by a single mode. While the inter-group anticrossings 
occasionally occur only between two neighboring groups for AB-stacked graphenes, they can 
be triggered among all groups for ABC-stacked graphenes. On the other hand, a signature 
characteristic of the ABC system is the appearance of the intra-group LL anticrossings in 
the energy region in which the LLs constituted by the quantized states on the sombrero¬ 
shaped subband surface densely appear. The low-lying magneto-electronic spectrum can 
be verified by STS and QBE experiments, while the confirmation of the entire spectrum 
can be achieved through optical spectroscopy measurements over a wide energy range, 
including absorption peaks, spectral widths and intensities of different optical transitions, 
as in Chapters 4.4, 5.5, and 6.5. 

7.2 Wave functions and density of states of Landau levels 

The density of states (DOS), reflecting the main characteristics of of the LL spectrum, is 
dehned as D{u)) = The i?o-dependent DOS exhibits many 

symmetric delta-function-like peaks attributed to the discrete LLs. The LL energies can 
be verihed by the energies of the sharp peaks in the differential conductance map of STS 
{dl/dV versus D), where the tunneling conductance (dl/dV) is proportional to the DOS at 
the tunneling electron energy [16611171] . For the DOS of the AA-stacked trilayer graphene, 
the energies of the zeroth LLs of the first, second and third groups are, respectively, eval¬ 
uated as —0.47, 0.01 and 0.46 eV (black, red and blue curves), as shown in Fig. 29(a). 
The peaks of every groups are located at energies followed by the relationship of monolayer 
graphene, i.e., cx A crossover of LLs in different groups leads to higher DOS 
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and stronger tunneling currents in experimental measurements. 

In AB-stacked trilayer graphene, the hrst group is followed by the square-root dependent 
sequence of monolayer graphene, and the second and third groups correspond to AB-stacked 
bilayer graphene, as shown in the Fig. 29(b). With half intensity, the peaks near the Fermi 
level are attributed to the symmetry-broken structure of AB-stacked trilayer graphene that 
leads to a considerable splitting (~10 meV) of the zeroth LLs in the hrst and second groups, 
while the spillting is hardly observable for the higher LLs under this held strength. On 
the contrary, the peak sequence of ABC-stacked trilayer graphene is neither like that of 
monolayer or bilayer graphenes, as shown in the Fig. 29(c). The hrst peak at the Fermi 
level, with intensity estimated to be three times than others, is a composition of three peaks 
originating from the Dirac points. Its intensity is proportionally increased with an increase 
of the number of graphene layers. On the other hand, the onset energies for the second 
and third groups are, respectively, evaluated as 0.34 eV and 0.40 eV, which approach to 
the crossover of the and S'g subbands under a weak magnetic held. Essential diherences 
of DOS are the spacings and heights of peaks, which can be verihed by STS. 

In addition, the real-space wave functions of LLs are observed by spectroscopic-imaging 
STM in 2D electron gas and topological insulator. From the measurements on the varia¬ 
tions of the local DOS in graphene planes, spectroscopic-imaging STM reveals the nodal 
structures of the wave functions attributed to the well-behaved LLs, where the zero-point 
numbers of the wave functions are identihed by htting the Fourier transformation (FT) of 
the LDOS data to the FT of the Landau wave functions. For the hybridized wave func¬ 
tions in the LL anticrossing regions, e.g., in AB- and ABC-stacked graphenes [^1731152] . 
spectroscopic-imaging STM provides a way to distinguish the main mode and the side 
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modes of the LLs. Also, it performs the verihcation on the LL splitting resulting from the 
symmetry broken down in AB-stacked graphene [TTHS^ . STS and STM directly measure 
the energy spectrum and the internal structure of wave functions, giving a useful iden- 
tihable picture of the Landau levels. The experimental results can be used as a tool to 
experimentally identify the stacking conhguration of a multilayer graphene. 

8 Concluding remarks 

The generalized tight-binding method is widely used to study the physical properties of 
massless and massive Dirac fermions in graphene-related systems. We review the magneto- 
electronic properties of multilayer graphenes, focusing on the energy spectra and the re¬ 
lated subenvelope functions in the AA, AB and ABC stacking conhgurations. Depending 
on the layer number, stacking conhguration and geometric symmetry, the LL spectra show 
feature-rich characteristics, including novel quantum numbers of the initial LLs, versatile 
held-dependent energies, high state degeneracy, semiconductor-metal transitions and inter- 
and intra-group anticrossings. Based on the sublattice dominance, the LLs of an A^-layer 
graphene can be clearly classihed into N groups, with the quantum numbers identihable by 
the derived dominant subenvelop functions. The degeneracy of LLs is broken by applying 
an external gate voltage to the graphene sheet. The onset energy and the held-dependent 
LL energies of each group are relevant to the zero-held band structure. The relationship 
among subenvelop functions of distinct sublattices is specihc for a LL under a given stack¬ 
ing conhguration. Via a perturbation theory, we can realize evolution of the subenvelop 
functions during the inter- and intra-group anticrossings. Few-layer systems have a rich 
variety of spectra with respect to diherent number of layers and stacking conhgurations. 
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Optical spectroscopy, STM and quantum Hall transport provide accuracy in measuring the 
low-lying magneto-electronic properties. The coupling of the inter- and intra-group anti¬ 
crossing LLs leads to unexpected physical properties, such as the extra absorption peaks 
and optical selection rules, and the variation of QHE conductivity plateaus and VHS in the 
DOS. In addition, spectroscopic-imaging STM can be used to distinguish the main mode 
and side modes of the hybridized LLs by measuring the variations of the LDOS. 

Magnetic quantization in AA-stacked multilayer graphenes results in the square-root 
LL energy dependence on the held strength and the quantum number, the monolayer char¬ 
acteristic, since the Dirac cones derived from the theoretical calculations and experiments, 
such as STS and ARPES, are preserved at the hexagonal corners, i.e., the K and K' points. 
For monolayer, the square-root dependence is limited to energies smaller than 0.4 70 , while 
with respect to each group in a multilayer case, the energy distribution is determined by 
whether the number of layers is odd or even. If the total number of layers N is odd, there 
is only one group situated close to the Fermi level, and other A^ — 1 groups are oppositely 
mirrored to each other with respect to the Fermi level. The Landau quantization mapped 
to the monolayer graphene is consistently presented once N is odd. In particular, the 
subenvelop functions with non-empty distribution have an identical relationship only for 
every other adjacent layer. In the even-A^ case, half the groups are above Ep={] and the 
other half lies below Ep = 0 with the mirror distribution described in accordance with 
the odd-A^ case; all layers are equivalent here. The quantum numbers for either case are 
identihable by the non-empty subenvelop functions of lattice A (B) of any layer. As a 
result, the zero-mode LLs correspond to the Dirac points and their energies are insensitive 
to the held strength. The absence of LL anticrossings in the energy spectra is consistent 
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with the fact that we can use a single mode to identify each LL due to the lack of coupling 
among LLs in the AA stacked conhguration. 

AB- and ABC-stacked graphenes both show a non-monotonic LL energy dependence 
and reveal a complex pattern of LL anticrossings due to the geometric symmetry. There is a 
total number of 4A^ Dirac-point LLs distributed over a narrow energy width (~10 meV) near 
Ep = 0. The LL spectrum of an A^-layer AB-stacked graphene consists of [(A^-|-l)/2] groups 
starting to form near Ep = 0 and another N — [(A^ -|-1)/2] groups away from Ep = 0. Once 
N is odd, one group of LLs especially exhibits a monolayer-like energy spectrum, in which 
the Landau states are contributed by electronic states localized on every other adjacent 
layer. The survival of monolayer properties is attributed to the Landau quantization of 
the massless Dirac quasiparticles. The degeneracy of LLs is broken in the odd-layer case, 
while it is protected in the even-layer case due to the inversion symmetry. Quantum 
transport experiments have verihed the variations in state degeneracies caused by the LL 
crossing and asymmetry-induced breaking through measurements of the Hall conductance 
plateaus. The dominating sublattices for energy subbands under magnetic quantization 
obey a regular rule. The two outermost layer sublattices are for the monolayer-like LLs, 
and the other different inner layer sublattices B{A) are for the other groups starting near 
Ep = 0 (away from Ep = 0). Furthermore, the anticrossing of LLs is attributed to the fact 
that the interlayer atomic interaction 73 couples the LLs in the 3 modulo, which prevents 
multi-mode LLs with certain identical modes from crossing each other in the Ho-dependent 
energy spectrum. 

For ABC-stacked graphenes, a versatility of magneto-electronic properties is observed 
when the number of layers is increased. With an onset energy near Ep = 0, there is always 
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one group of LLs responsible for the magnetic quantization of the surface localized states. 
The phenomenon of LL anticrossings is strikingly pronounced in ABC-stacked graphenes, 
in contrast to the few anticrossings observed in AB-stacked graphenes and the absence of 
anticrossings in AA-stacked graphenes. In the energy spectrum, the LL evolution under 
a magnetic held reveals a complex pattern of two kinds of anticrossings by intergroup 
and intragroup LLs, and some peculiar LLs attributed to the Landau quantization of the 
sombrero-shaped energy subbands. The intergroup anticrossings appear for any two groups 
under a sufficiently strong magnetic held. On the other hand, intragroup LLs frequently 
anti-cross each other in a small energy region where the peculiar LLs, with energies inversely 
dependent on i?o-held strength, are arranged in the order of decreasing quantum numbers 
when the energy is increased. However, the unusually sequenced LLs and the intragroup 
anticrossings start to disappear at a critical held strength, in which Be is higher for a 
system with more layers. In addition to the nonvertical hopping integral (3^ between two 
adjacent layers, the two hopping integrals (32 and (3^ between two next-neighboring layers 
also take the same role and combine the LLs in the 3 modulo in ABC-stacked graphenes. 
This is attributed to the specihe relationship between the subenvelop functions and the 
interlayer atomic interactions. 

A generalization of magneto-electronic properties is valuable in understanding other 
physical properties, such as the mechanisms behind magneto-optical excitations. The in¬ 
tragroup transitions between well-dehned LLs satisfy the particular selection rule An = ±1 
in the AA, AB and ABC stacking conhgurations, whereas the rule for intergroup transitions 
depends on the conhgurations. The former is also applicable for the LLs of Dirac quasi¬ 
particles in the Graphene-like 2D materials, e.g., M 0 S 2 [11411115] and Silicene [11411116] and 
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the topological insulator |11811121j . The intergroup transitions are forbidden due to the 
derived zero electric dipole moment in the AA-stacked system. Also, for the AB-stacked 
graphenes with odd number of stacked layers, excitations from the quantized LLs of the 
linear bands to those of the parabolic bands are forbidden. However, considering intergroup 
transitions for the other cases, the optical selection rules are defined by the concept of dipole 
transition. Furthermore, the hybridization of the LLs induces new optical selection rules for 
modulo 3 in addition to the An = ±1 applied to transitions between the well-defined LLs. 
Moreover, the rich magnetic quantization leads to the complicated electronic excitations, 
mainly owing to the strong competition between the transverse Lorentz motion and the 
longitudinal Coulomb oscillation. There are a lot of low-frequency magneto-plasmons with 
a novel momentum-frequency relationship. The AA- and AB-stacked bilayer graphenes, 
respectively, exhibit the 2D-like plasmon and the discrete inter-LL plasmons. For mul¬ 
tilayer graphene, the main characteristics of magneto-plasmons, existence, intensity and 
frequency, will be significantly affected by the stacking configuration and the layer number. 

Ever since graphene has first been successfully fabricated, its unique physical prop¬ 
erties have drawn intensive attentions of scientists that are devoted to researching new 
two-dimensional materials for possible applications in next-generation nanoelectronic de¬ 
vices and fundamental scientific research. Recently, 2D materials like M 0 S 2 [22211223] and 
silicene |224fl226] . arranged within a buckled honeycomb lattice, have been synthesized and 
present an intrinsic bandgap |222[|223]I225[I226] . The electronic dispersion resembles that 
of relativistic Dirac fermions. Likewise, the generalized tight-binding model can cope with 
the sp^ hybridizations in the buckled structure and handle the influence of external helds; 
it can be adapted to different geometric structures, including bounded, folded, curved, and 
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stacked configurations. 
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Figure Captions 

Figure 1. (a) Honeycomb lattice structure of monolayer graphene. Sublattices A and 
B are shown by black and red colors, respectively. The primitive unit cell is depicted by 
a gray diamond, where a indicates the hopping integral between two sublattices A and B 
and b' represents the C-C bond length, (b). The Brillouin zone of honeycomb lattice and 
some highly symmetric points, (c). The magnetically induced rectangle unit cell under a 
uniform magnetic held B = Bqz, perpendicular to the graphene plane. 

Figure 2. Geometric structure of AA-stacked graphene. The interlayer atomic interac¬ 
tions are illustrated in the right panel. 

Figure 3. Geometric structure of AB-stacked graphene. The interlayer atomic interac¬ 
tions are illustrated in the right panel. 

Figure 4. Geometric structure of ABG-stacked graphene. The interlayer atomic inter¬ 
actions are illustrated in the right panel. 

Figure 5. Band structure of monolayer graphene along F —)■ iF —)■ M direction indicated 
in Fig. 1(b). 

Figure 6. (a). Landau levels of monolayer graphene at Bo = 20 T. The quantum 

numbers = 0,1, 2, 3... are counted from the Fermi level, Ep = 0, to higher conduction 
and valence energies, (b) and (c). The envelope functions of two sublattices A and B are 
plotted with respect to the Landau levels shown in (a). 

Figure 7. Energy dependence of Landau level on the square root of (a) magnetic held 
strength and (b) quantum number. 

Figure 8. Band structure of monolayer graphene under a modulated magnetic held with 
the modulation held strength Bm = 20 T and the period length Rm = 500. The center of 
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the oscillating Landau subbands is localized at ky = 2/3. Red lines indicate the energies of 
the Landau levels under a uniform field of 20 T. 

Figure 9. (a) Quasi-Landau levels of flat graphene ribbon. Oscillating Landau subbands 
of (b) curved graphene ribbon and (c) carbon nanotube. The parameters W, 6 and R 
represent the width of ribbon, the arc angle of curved ribbon and the radius of tube, 
respectively. 

Figure 10. Zero-field band structure of AA-stacked graphenes for (a) odd number of 
layers: monolayer and trilayer, and (b) even number of layers: bilayer and tetralayer. 

Figure 11. (a) Landau levels of AA-stacked trilayer graphene under Bq = 20 T. For 
each level, the envelop functions of sublattices A^, A^, B^, A^ and B^ for the 1/6 

localization state are shown in (b), (c), (d), (e), (f) and (g), respectively. 

Figure 12. Landau-level spectrum of AA-stacked trilayer graphene. The Fermi level is 
plotted as a bold wiggling curve, an oscillation around the zero energy as a function of the 
magnetic field strength. 

Figure 13. Landau-level spectrum of AA-stacked bilayer graphene. The Fermi level is 
plotted as a bold wriggling curve, an oscillation around the zero energy as a function of 
the magnetic held strength. 

Figure 14. Landau-level spectrum of AA-stacked tetralayer graphene. The Fermi level 
is plotted as a bold wriggling curve, an oscillation around the zero energy as a function of 
the magnetic held strength. 

Figure 15. Metal-semiconductor transitions during a variation of the magnetic held for 
AA-stacked (a) bilayer and (b) trilayer graphenes with/without an applied gate voltage. 

Figure 16. Band structure of AB-stacked (a) bilayer and (b) trilayer graphens in the 
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absence of external fields. The inserts of (a) and (b) show the zoomed-in view near the 
Fermi level. 

Figure 17. (a) Landau levels of AB-stacked trilayer graphene under Bq = 20 T. For 
each level, the envelop functions of sublattices A^, A^, B^, and B^ for the 1/6 

localization state are shown in (b), (c), (d), (e), (f) and (g); those for 2/6 localization state 
are shown in (h), (i), (j), (k), (1) and (m), respectively. 

Figure 18. Landau-level spectrum of AB-stacked trilayer graphene. 

Figure 19. Landau-level spectrum of AB-stacked bilayer graphene. 

Figure 20. Landau-level spectrum of AB-stacked tetralayer graphene. 

Figure 21. Band structure of ABC-stacked (a) trilayer and (b) tetralayer graphenes in 
the absence of external helds. The inserts of (a) and (b) show the zoomed-in view near the 
Fermi level, where the crossing points are the Dirac points. 

Figure 22. (a) Landau levels of ABC-stacked trilayer graphene under Bo = 20 T. For 
each level, the odd-indexed envelop functions of sub lattices Al, i?/, A^, B^, A^ and B^ for 
the 2/6 localization state are shown in (b), (c), (d), (e), (f) and (g), respectively. 

Figure 23. Landau-level spectrum of ABC-stacked trilayer graphene. 

Figure 24. (a) Landau-level anticrossing pattern of ABC-stacked trilayer graphene. For 
the states marked by circles in (a), their envelop functions of sublattices A^, B^, Al, B^, 
Al and B^ for the 2/6 localization state are plotted in (b), (c), (d), (e), (f), and (g), 
respectively. 

Figure 25. Landau-level spectrum of ABC-stacked tetralayer graphene. 

Figure 26. Landau-level spectrum of ABC-stacked pentalayer graphene. The main 
characteristics of the hrst to hfth group are plotted in (a), (b), (c), and (d). 
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Figure 27. Energies of Dirac-point related Landau levels of ABC-stacked trilayer graphene 
are plotted for i?o=25 T based on their dependence on the number of stacked layers. 

Figure 28. (a) Landau levels of the second group of ABC-stacked trilayer graphene are 
plotted based on the full tight-binding model (black) and the minimum model (red). Red 
cross symbols indicate the highest Landau levels which are fully occupied under the critical 
helds Bg. (b) Linear relationship between Bg and Ug is plotted, where Ug appears with a 
constant period of l/i?o- 

Figure 29. Density of states of Landau levels for (a) AA-, (b) AB-, and (c) ABC-stacked 
trilayer graphenes. 
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Figures 10 and 11 
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Figure 12 
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Figure 13 
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Figure 15 
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Figures 16 and 17 
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Figure 18 
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Figures 21 and 22 
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Figure 23 


96 













A 

It 

JL 

JL 

4- 

Jl 

Jl 

jii_ 

JL 

jL 

Hi 

(c) 

b: 

TIL 

Tit 

TIL 

jL 

JiV- 

jvL 

AIL 

jiL 

jii 

-AH— 

(d) 

a: 

_/VV\_ 

_/vv>^ 

_/W^ 

-TWL_ 

/Wv 

J»L 

JL 

JL 

JL 

/A 

(e) 

1 

x 

i: 

X 

± 

Jl. 

-A_ 

A 

-A- 

/H 

(0 

a: 

— _ 

—WV— 

— 

— 

_ 

-JVV— 

JL 

jA 

/A 

(g) 

b: 

JW\_ 

TVA 

-TY^ 

JL 

JL 

JL 

JL 

A 

A 

iL 

4 

It 

fA 

/i 

A 











10 20 30 40 50 0.31 1/3 


1/3 1/3 


1/3 1/3 


1/3 


B(T) Ai,(BD/2RB 


Figure 24 
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Figure 25 
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Figure 26 
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Figures 27 and 28 
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Figure 29 
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